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ABSTRACT 

We present our X-ray imaging spectroscopic analysis of data from deep Suzaku and XMM-Newton 
Observatory exposures of the Virgo Cluster elliptical galaxy NGC 4649 (M60), focusing on the abun- 
dance pattern in the hot interstellar medium (ISM). All measured elements show a radial decline in 
abundance, with the possible exception of Oxygen. We construct steady state solutions to the chemi- 
cal evolution equations that include infall in addition to stellar mass return and SNIa enrichment, and 
consider recently published SNIa yields. By adjusting a single model parameter to obtain a match to 
the global abundance pattern in NGC 4649 we infer that introduction of subsolar metallicity external 
gas has reduced the overall ISM metallicity and diluted the effectiveness of SNIa to skew the pattern 
towards low a/Fe ratios, and estimate the combination of SNIa rate and level of dilution. Evidently, 
newly-introduced gas is heated as it is integrated into, and interacts with, the hot gas that is already 
present. These results indicate a complex flow and enrichment history for NGC 4649, reflecting the 
continual evolution of elliptical galaxies beyond the formation epoch. The heating and circulation of 
accreted gas may help reconcile this dynamic history with the mostly passive evolution of elliptical 
stellar populations. In an appendix we examine the effects of the recent updated atomic database 
AtomDB in spectral fitting of thermal plasmas with hot ISM temperatures in the elliptical galaxy 
range. 

Subject headinqs: galaxies: abundances, galaxies: elliptical and lenticular, galaxies: individual (NGC 
4649), galaxies: ISM 


1. INTRODUCTION 

There are two basic approaches to studying the forma- 
tion and evolution of elliptical galaxies. One may directly 
examine the assembly history of the baryonic compo- 
nent by observing how the space density and morpholog- 
ical demographics of the population of elliptical galaxies, 
and their progenitors, develop over time. Alternatively 
one may conduct detailed investigation of those spectro- 
photometric properties that reflect their histories, and 
the scaling relations among these properties. While sur- 
veys at a range of redshifts indicate significant growth in 
both the number and sizes of ellipticals, “archeological” 
investigation finds that elliptical galaxy stellar popula- 
tions are mostly in place at high redshifts (z > 2 for 
massive systems) and passively evolve thereafter. Reso- 
lution of this apparent paradox is crucial for validation 
and elucidation of the prevailing hierarchical assembly 
paradigm of galaxy formation, where dynamical evolu- 
tion and the triggering of star formation are presumed 
to be interconnected. 

The archeological approach traditionally relies on op- 
tical photometry and spectroscopy of the stellar popu- 
lation; however, the hot interstellar medium (ISM) pro- 
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vides a rich complementary site of diagnostic data that 
is accessible by means of X-ray observation (Loewenstein 
& Davis 2010, hereafter Paper I, and references therein; 
Pipino & Mattecucci 2011, hereafter PM11). The phys- 
ical properties of the ISM in giant elliptical galaxies re- 
flect the distinctive history and nature of these systems 
and, as such, markedly differ from those in the ISM of 
spiral galaxies such as the Milky Way. The most strik- 
ing contrast is that, while mass loss from evolved stars 
is a primary source of gas in both spirals and ellipticals, 
most of the mass return in the latter is promptly heated 
to the high temperatures corresponding to the stellar ve- 
locity dispersion and does not currently participate in an 
ongoing star-gas cycle (Mathews & Brighenti 2003). 

This distinct ecology provides a repository of informa- 
tion about processes in the distant and recent past of 
elliptical galaxies. The ISM is more responsive to ener- 
getic events than the stellar population which is, to first 
order, passively evolving. Evidence of feedback processes 
subsequent to the establishment of the stars and the nu- 
clear supermassive black hole (SMBH) can only be found 
in the ISM. The ISM mass within elliptical galaxies is a 
small fraction of the stellar mass return integrated over 
the post-star-formation history of an elliptical galaxy, 
implying the existence of some steady and/or episodic 
means of gas removal into the intergalactic, or some cir- 
cumgalactic, medium. A successful gasdynamical model 
for ellipticals must explain ISM that generally are, nev- 
ertheless, sufficiently massive to rule out persistent su- 
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personic or transonic galactic winds, and also display a 
large scatter in gas-to-stellar mass ratio (M gas /M sta rs ; 
Mathews et al. 2006 and references therein). While ram 
pressure stripping may be important at times for some el- 
lipticals in rich clusters, large scale flows driven by Type 
la supernovae (SNIa) and active galactic nuclei (AGN) 
likely dominate the redistribution of the ISM and its re- 
moval from the galaxy. 

The amount of energy associated with SNIa exploding 
at the estimated rate in ellipticals is sufficient to drive a 
galactic wind. Although effects such as depletion through 
the formation of dust (PM11) may reduce the discrep- 
ancy, the high expected ISM Fe abundance is at odds 
with X-ray observations (e.g., Paper I, PM11), and a 
certain amount of fine-tuning in the SNIa rate as a func- 
tion of time is required for SNIa-driven winds to explain 
the Aigas/Mstars scatter without driving out virtually all 
of the gas and rendering ellipticals undetectable as dif- 
fuse soft X-ray sources. Feedback from active galactic 
nuclei (AGN), implicated in quenching star formation in 
ellipticals (Schawinski et al. 2007; Cattaneo et al. 2009) 
and establishing the scaling relations between SMBH 
mass and stellar mass or velocity dispersion (e.g., De- 
Buhr, Quataert, & Ma 2011 and references therein), may 
also drive galactic flows at later times. AGN interaction 
with hot ISM is evident in Chandra X-ray Observatory 
observations of several ellipticals (McNamara & Nulsen 
2007; Nulsen et al. 2009). AGN feedback is fundamen- 
tally self-regulating and intermittent: the SMBH is fu- 
eled by an initial inflow that is subsequently reversed 
as the AGN powers up, thus cutting off the supply of 
gas and enabling the cycle to restart as gas once again 
flows inward towards the now-dormant SMBH. There 
has been great progress in implementing and applying 
hydrodynamical simulations that include various feed- 
back prescriptions (Mathews, Brighenti, & Buote 2004; 
Brighenti et al. 2009; Ciotti, Ostriker, & Proga 2010), 
however it remains to be seen whether the full range 
of ISM observables - X-ray luminosities, tempera- 
tures, metallicities and abundance patterns - and their 
galaxy-to-galaxy variations can be accurately and self- 
consistently modeled without additional, and perhaps 
fundamental, adjustment (Mathews & Brighenti 2003). 
In addition, many ellipticals are embedded in an ex- 
tended hot intergalactic or circumgalactic medium, and 
may interact with their environment in a variety of ways 
(Mulchaey & Jeltema 2010). In particular ellipticals 
may accrete some combination of primordial gas and gas 
ejected from winds at earlier epochs (Pipino & Matteucci 
2004; Dave et al. 2011). Such an external medium may 
also confine outflows. 

Recent X-ray measurements of elemental abundances 
in the hot ISM of ellipticals beg consideration in deter- 
mining the future direction of these models, and con- 
straining more general theories of the chemical evolution 
of these systems. The ISM abundance pattern may be 
a particularly sensitive diagnostic of dynamical processes 
in ellipticals. If the sole source of interstellar gas is stellar 
mass loss, the ISM abundance pattern in ellipticals will 
reflect that in evolved stars found to have [a/Fe] s t ars > 0 
(\a/Fe\ stars is defined as the log of the abundance ra- 
tio, with respect to Fe, of elements primarily produced 
through a capture to Fe - relative to the solar ratio). 
Because the internal injection and flow of energy, and 


the various mechanisms of mass exchange with the ex- 
ternal environment, each imprint distinctive departures 
from the baseline ISM abundance pattern determined by 
local stellar mass loss, this pattern may be analyzed to 
probe for signatures of these processes. 

Measurement of elliptical galaxy abundance patterns 
encompassing a broad range of elements, and extending 
to large radii may be made with the Suzaku Observatory, 
facilitated by the low internal background and relatively 
sharp energy resolution of the XIS CCD detectors. In 
Paper I we derived the abundance pattern in the ellip- 
tical galaxy NGC 4472 from analysis of Suzaku spectra, 
supported by analysis of co-spatial XMM-Newton Obser- 
vatory spectra. Application of simple chemical evolution 
models to these data, led us to conclude that the abun- 
dances may be explained by a combination of a-element 
enhanced stellar mass loss and direct injection of ejecta 
from SNIa exploding at a rate ~ 4 — 6 times lower than 
the standard value. In addition, we discovered abun- 
dance anomalies in the sense that no published set of 
SNIa yields could simultaneously reproduce the inferred 
Ca and Ar, and Ni abundances; and (confirming prior re- 
sults summarized in Paper I) in the sense that standard 
core collapse nucleosynthesis models evidently overpro- 
duce O by ~ 2. 

In this paper we adopt a broadly similar approach in 
our investigation of NGC 4649 (M60). As is the case for 
NGC 4472, NGC 4649 is a giant elliptical galaxy in the 
Virgo cluster with an old stellar population enhanced in 
a-elements - but with several notable differences both 
optically and in X-rays. In particular NGC 4649 has a 
substantial major axis rotation (Brighenti et al. 2009, 
and references therein), and is considerably more com- 
pact in X-rays. We have adjusted our spectral analysis 
procedures in response to the distinctive X-ray charac- 
teristics of NGC 4649, and have revised and expanded 
our models to consider possible episodes of inflow, and a 
recently published set of SNIa yields. 

We detail our data reduction and spectral analysis pro- 
cedures in Sections 2 and 3, where we present our derived 
NGC 4649 hot ISM thermal and chemical properties and 
their radial variation. In Section 4 we focus on interpret- 
ing the global abundance pattern in the context of the 
relative contributions of metal enrichment from stellar 
mass return, SNIa, and inflow of extragalactic material 
using steady state solutions to the equations of chemical 
evolution. Section 5 includes a summary of our conclu- 
sions, and discusses possible implications of our results. 
In an appendix we examine the effects of the recent up- 
dated atomic database AtomDB in spectral fitting of el- 
liptical galaxy hot ISM such as NGC 4649. 

2. CONSTRUCTION OF NGC 4649 SPECTRA AND 
ASSOCIATED FILES 

2.1. Suzaku Spectral Extraction and Preparation 

NGC 4649 was observed with Suzaku (Mitsuda et al. 
2007) between 2006-12-29 and 2007-01-04 (ObsID 
801065010), at which time three co-aligned, 17.8' x 17.8' 
held-of-view X-ray Imaging Spectrometer (XIS) CCD 
cameras (Koyama et al. 2007) - two front-illuminated 
(FI: XIS0 and XIS3) and one back-illuminated (BI: 
XIS1) - were operational, each XIS in the focal 
plane of an X-ray Telescope (XRT) with a 2' half- 
power diameter (Serlemitsos et al. 2007). Observations 
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were conducted utilizing the space-row charge injection 
(SCI) technique that reverses the degradation in en- 
ergy resolution caused by accumulated radiation dam- 
age (Nakajima et al. 2008). We initiate our data re- 
duction with the unfiltered event files. These data un- 
derwent Version 2.0.6.13 pipeline processing on 2007- 
OS- 17 that enables one to properly account for the ef- 
fect of SCI on instrument characteristics and perfor- 
mance (Uchiyama et al. 2008). We reprocess the unfil- 
tered event files by hand in order to apply updated cal- 
ibration data and software, following the procedures de- 
scribed in Paper I. We recalculate PI values and grades, 
select event grades (0, 2, 3, 4, 6) that correspond to X- 
ray photon events, filter on pixel status (eliminating bad 
charge transfer efficiency columns, and rows invalidated 
by the charge injection process) and select good time 
intervals (GTI) based on pointing, data and telemetry 
rates, SAA proximity (“SAA.HXD = 0,T_SAA_HXD > 
436”), and proximity to the earth’s limb and illumi- 
nated Earth (“ELV > 5, DYE.ELV > 20”). In addi- 
tion, telemetry-saturated frames and calibration source 
photons are screened out; and, hot and flickering pixels 
are removed, and we accept only GTI where the revised 
geomagnetic cut-off rigidity COR2 > 4, thus eliminat- 
ing intervals with the highest particle background level 
(Tawa et al. 2008) without compromising overall statis- 
tical accuracy. 5x5 event files are converted to 3 x 3 
mode format, and merged with the 3x3 event files. The 
exposure times in the cleaned event files are 216.6 (216.5) 
ks for the XIS0 and XIS3 (XIS1) spectra. 

We extract spectra from the inner 8' (~ 6.5 R e = 36.2 
kpc at the NGC 4649 distance of 15.6 Mpc; Tonry et 
al. 2001, where R e is the optical half-light radius from 
Bregman, Temi, & Bregman 2006) in 2' and 4' circu- 
lar annuli, centered on the NGC 4649 optical nucleus 
(a = 12/i43m40.0s, <5 = +11°33'10") that very closely 
corresponds to the X-ray peak in the Suzaku image. 
The spectral redistribution matrix files (rmf) are gen- 
erated using xisrmfgen version 2009-02-28. The rmf 
and spectral files are binned to 2048 channels. The ef- 
fective area function files (arf) for the source spectra 
are generated by the xissimarf gen version 2009-01-08 
Monte Carlo ray-tracing program (Ishisaki et al. 2007) 
with 400000 simulation photons per energy bin, and an 
input source fits image file generated from a /3-model fit 
to the background-subtracted X-ray surface brightness 
profile extracted from archival Chandra data (/3 = 0.467, 
core radius 3.5"). Spectra from the FI chips, XIS0 and 
XIS3, are co-added and a weighted XISO+3 response 
function calculated from their respective rmf and arf 
files. Source spectra are grouped into bins with a min- 
imum of 15 cts, and we derive best-fit parameters and 
confidence levels using y 2 statistics. Errors correspond 
to 90% confidence limits. 

The XIS background consists of Non- X-ray (instru- 
mental and charged particle) Background (NXB), Galac- 
tic X-ray Background (GXB), and (extragalactic) Cos- 
mic X-ray Background (CXB) components. Since NGC 
4649 fills the Suzaku field of view, we estimate and sub- 
tract the NXB component and include additional com- 
ponents in our spectral fits to account for the GXB and 
CXB. The NXB component is estimated from observa- 
tions of the night earth taken in SCI mode within 150 
days of the starting or ending dates of our observation 


Table 1 

Suzaku Count Rates 


region 

XISO+3 

XIS1 

0-2' 

0.399 (0.008) 

0.299 (0.011) 

2 - 4' 

0.154 (0.058) 

0.127 (0.069) 

4-6' 

0.0840 (0.18) 

0.0676 (0.21) 

6-8' 

0.0583 (0.18) 

0.0522 (0.22) 

4-8' 

0.151 (0.22) 

0.124 (0.25) 

2-6' 

0.232 (0.10) 

0.189 (0.125) 


Note. — Total counts rates in cts s — 1 in the energy band used for 
spectral fits (0.4-4. 0 keV for the 6 — 8' annulus, 0.4-7. 0 keV for all 
others). Values in parentheses represent the fractions of total counts 
in the NXB. 

using xisnxbgen version 2008-03-08. The NXB event 
list in that time interval undergoes the identical screen- 
ing as the source data, is sorted by geomagnetic cut-off 
rigidity, and weighted according to the cut-off rigidity 
distribution in the source event file (Tawa et al. 2008). 
The estimated NXB spectra include only those events 
collected in the regions on the detector from which the 
source spectra are extracted. Since the CXB and GXB 
are spatially distributed distinctly from that of the NGC 
4649 X-ray emission, a separate arf file is generated using 
2000000 simulation photons per energy bin from a uni- 
form source of radius 20' and applied to the background 
in spectral fits. Combined XISO+3 NXB spectra and re- 
sponse functions are constructed as described above. 

We display the NXB-subtracted source count rates in 
each annulus, and the fraction of the total counts in the 
NXB for the corresponding bandpass used in spectral 
fitting (see below), in Table 1 . 

2.2. XMM-Newton EPIC Spectral Extraction and 
Preparation 

NGC 4649 was observed with XMM-Newton 
(Jansen et al. 2001) for ~54 ks on 2001-01-02 (Ob- 
sID: 0021540201), and for ~92 ks on 2007-12-19 (ObsID: 
0502160101). We downloaded the data from the High 
Energy Science Archive Research Center (HEASARC) 
data archive 3 and processed them with Science Analysis 
System (SAS) version 11.0.1, utilizing the techniques de- 
scribed in Kuntz & Snowden (2008) and Snowden et al. 
(2008) through application of the XMM-ESAS (XMM- 
Newton Extended Source Analysis Software) suite of 
procedures. EPIC-MOS and -PN CCD camera photon 
event files were generated in the standard way with the 
tasks emchain and epchain, respectively 4 . Files were 
processed with FLAG == 0, and screened to retain 
only those events with PATTERN <= 12 for the MOS 
and PATTERN <= 4 for the PN. 

The event files for each observation was screened 
for flares by examining their lightcurves using the 
xmmlight_clean scripts (Bauer et al. 2008). The use- 
ful exposure times, after cleaning, are ~49 ks (MOS) 
and ~41 ks (PN) for the ObsID 0021540201 event lists, 
and ~73 ks (MOS) and ~70 ks (PN) for the ObsID 
0502160101 event lists. XMM-Newton observations may 
be significantly contaminated by soft protons originat- 

3 http:/ /heasarc. gsfc.nasa.gov/docs/archive.html 

4 http:/ /xmm.esac.esa.int/external/xmm_user_support/ 
documentation/sas_usg/USG/ 
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Table 2 

XMM-Newton EPIC Count Rates 


Table 3 

XMM-Newton RGS Count Rates 


region 

MOS1 

Exposure MOS2 Exposure 

PN 

Exposure 

0 

- 2' 

0.518 

45.21 

ObsID 0021540201 
0.507 45.54 

1.797 

49.28 

2 

- 4' 

0.173 


0.185 

0.672 


4 

-6' 

0.103 


0.103 

0.192 


6 

-8' 

0.067 


0.059 

0.313 


0 

- 2' 

0.504 

67.24 

ObsID 0502160101 
0.482 66.46 

1.782 

58.25 

2 

- 4' 

0.154 


0.142 

0.447 


4 

-6' 

0.103 


0.080 

0.201 


6 

-8' 

0.056 


0.057 

0.128 



Note. — Total counts rates in cts s — 1 in the energy band used for 
spectral fits (0.4-5. 5 keV for the PN and 0.3-5. 5 keV for the MOS). The 
exposure time is in units of ks. 


ing in the solar wind (Kuntz & Snowden 2008). Al- 
though the lightcurve filtering process largely removes 
soft proton flares, we determine the level of residual soft 
proton contamination using the Fin_over_Fout method 
(Molendi, De Luca, & Leccardi 2004). We find the level 
of contamination to be minor for the MOS detectors, and 
negligible for the PN detector, during the first (2001) 
observation. However for the second (2007) observation, 
although the two MOS cameras again experienced only 
minor contamination, the PN camera was significantly 
affected. 

The extracted spectra are binned so that each channel 
has a minimum of 50 counts, and we consider an energy 
range of 0.3-5. 5 keV for the MOS detectors and 0.4-5. 5 
keV for the PN detector. We extracted spectra from con- 
centric annuli with identical centers and (2') widths to 
those of the Suzaku extraction regions. The spectral ex- 
traction was performed using the XMM-ESAS tasks mos- 
spectra and pn-spectra for the relevant detectors. These 
scripts extract the spectra, and generate the instrument 
response (rmf and arf) matricies. The tasks mos-back 
and pn-back are used to model the instrumental and par- 
ticle backgrounds (NXB) for each of the extracted spec- 
tra. 

The source count rates in each annulus are displayed 
in Table 2. 

We use the the HEASARC background tool 5 to extract 
the ROSAT All-Sky Survey (RASS) spectrum from an 
annular region, with an inner radius of 1° and an outer 
radius of 2°, centered on the X-ray peak of NGC 4649 - 
assuming that the average in the annulus is a fair rep- 
resentation of the CXB and GXB at the position of the 
galaxy. This data is simultaneously fit with the EPIC 
data. 


2.3. XMM-Newton RGS Spectral Extraction 

We extract the RGS spectra for this galaxy using the 
standard SAS procedures, applying the rgsproc task to 
generate filtered event files and spectra for both observa- 
tions (OBSID 0021540201 and 0502160101). We exclude 
some short time intervals where we detected flares in the 
background light curves, resulting in the exposure times 
listed in Table 3. 

We use the standard background products and gener- 

5 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl 


region 

RGS1 

Exposure 

RGS 2 

Exposure 

Order 




ObsID 0021540201 



0 

- R5 

0.064 

52.56 

0.100 

51.08 

i 

0 

- 1?5 

0.026 


0.028 


2 




ObsID 0502160101 



0 

- 1.'5 

0.051 

74.73 

0.090 

74.78 

1 

0 

- 1?5 

0.026 


0.027 


2 


Note. — Total counts rates in cts s 1 for the RGS spectra for each 
observation. The exposure time is in units of ks. 


ate RGS response matrices using the task rgsrmfgen that 
are technically applicable to spectra of point sources in 
the RGS held of view. To account for the mixing of spa- 
tial and spectral information from extended sources, the 
matrices are convolved with an image of NGC 4649 (with 
point-sources removed) created from Chandra ACIS-S 
data (OBSID: 8182) using the FTQQL rgsrmfsmooth. The 
image is adaptively smoothed with minimum of 75 counts 
to reduce the noise prior to the response matrix convo- 
lutions. 


3. SPECTRAL ANALYSIS 
3.1. Suzaku Spectral Analysis 

We employ Xspec version 12.7 to simultaneously fit 
the Suzaku XISO+3 and XIS1 NGC 4649 spectra. Our 
baseline source model consists of an absorbed single- 
temperature thermal plasma, plus a 7.5 keV thermal 
bremsstrahlung component to account for the LMXB 
(Irwin et al. 2003; Humphrey & Buote 2008). Temper- 
atures and heavy elemental abundances (C, N, O, Ne, 
Mg, Al, Si, S, Ar, Ca, Fe, Ni) in the former com- 
ponent are estimated using the vapec model that, in 
this latest version of Xspec, incorporates the updated 
AtomDB v2.0.1 atomic database. We note that, at 
elliptical galaxy X-ray temperatures, spectral fits are 
significantly impacted by changes with respect to the 
AtomDB vl.3.2 database used in earlier versions of 
Xspec (see Appendix A). Therefore, any detailed com- 
parison with previous results (Randall, Sarazin, & Irwin 
2004, 2006; Tozuka & Fukazawa 2008; Humphrey et al. 
2008; Diehl & Statler 2008; Nagino & Matsushita 2009) 
are of limited utility, although we note a general < 4' 
ISM temperature profile consistency if we use AtomDB 
1.3.2. Addition of a second vapec model (with identi- 
cal abundances) does not improve the fits, and is not 
further considered. We fix the redshifts of the source 
components at the value corresponding to the NGC 4649 
distance. We apply the Wilms, Allen, & McCray (2000) 
Tuebingen-Boulder ISM absorption model (tbabs), and 
values in excess of the Galactic value of 2.2 x 10 20 cm' 2 
(Dickey & Lockman 1990) are considered as described 
below. We adopt the Asplund et al. (2009) solar abun- 
dance standard, newly implemented in Xspec 12.6, fix- 
ing C and N at their solar values. Since we find that only 
O, Ne, Mg, Si, S, and Fe are (formally) well-constrained, 
Al:Mg, Ca:Ar:S, and Ni:Fe are fixed at their solar ratios 
(Murakami et al. 2011). The models include a constant 
XIS1/XIS03 multiplicative factor to account for flux cal- 
ibration offsets that, when fitted for, ranges from 1.0- 
1.08. 
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We include four astrophysical background components 
for all extraction regions, although these are not always 
formally required for the innermost annuli where the 
contribution of one or more of these is small or negli- 
gible. We model the CXB by an index 1.4 power-law 
(Kushino et al. 2002) , and the GXB by the sum of solar- 
abundance thermal plasma (apec) models representing 
the Milky Way halo (MWH) and Local Hot Bubble 
(LHB) (Hayashi et al. 2009; Nagino & Matsushita 2010). 
The CXB and MWH are absorbed by the Galactic col- 
umn, the LHB is unabsorbed. We also include an ad- 
ditional extended, 0.5 solar abundance apec component 
(Loewenstein, Valinia, & Mushotzky 2001) to model the 
Virgo ICM (see, also Buote & Fabian 1998, Randall et al. 
2006, 2008, Ji et al. 2009 - hereafter Ji09, and Bogdan 
& Gilfanov 2011, for evidence of emission from a second, 
hotter medium). We allow normalizations for these com- 
ponents in each annuli to freely vary in general; however, 
the MWH/LHB flux ratio is assumed to be constant (see 
below). We tie the background XISl-to-XISO+3 ratio 
to the value in the source model. Separate arfs are ap- 
plied to background and source model components as 
described above. 

Because the GXB is most prominent in the outermost 
annuli, we fit the 6 — 8' annulus (simultaneously with the 
RASS spectrum) first to determine the LHB, MWH, and 
ICM temperatures. The NXB dominates above 4 keV in 
this annulus, so we fit over the 0.4-4 keV bandpass (all 
other spectra are fit over 0.4-7 keV). The column den- 
sity is fixed at the Galactic value, however post facto fits 
with the higher absorption determined from the inner- 
most annulus (see below) yield virtually identical results. 
The background-component temperatures in the spectra 
extracted from other regions are fixed at these outer re- 
gion best-fit values: fcT(LHB) = 0.105, fcT(MWH) = 
0.187, fcT(ICM) = 1.69, as is the MWH/LHB flux ratio 
(Murakami et al. 2011). Relaxing these constraints did 
not significantly impact the best-fit ISM parameters. 

We display the minimized (reduced) y 2 values, and 
source model temperatures abundances, and column den- 
sities, for each annuli in Table 4; and plot spectra, best-fit 
models (and their components), and data residuals to the 
best-fit model for the 2! annular partitioning of the 8' re- 
gion in Figures 1-4. We note the following details. We 
find degeneracies in some fits among the normalizations 
of the CXB, ICM, and LMXB hard component normal- 
izations (thus the latter component does not appear in 
Figures 1 and 3); however, these do not introduce sig- 
nificant errors in the hot ISM abundance or tempera- 
ture determinations. Best-fit model parameters of inter- 
est are insensitive to whether thermal bremsstrahlung or 
power-law models are adopted for LMXB components. 
We adopt the former since in this case fits are somewhat 
less sensitive to assumptions about the absorption. In 
the 0 — 2' region, we found systematic positive residu- 
als at ~ 1.23-1.24 keV that can be well fit by a nar- 
row ( i.e . unresolved by Suzaku ) Gaussian emission line. 
We include this component in the fits plotted and re- 
ported here. This excess may be a sign of lingering in- 
adequacies in calibration (via the instrumental response 
or background) or in the plasma code; and, modeling 
it only effects the quality of fit (a reduction of y 2 by 
~ 100) and does not alter the best-fit parameters. Ex- 
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Figure 1. 6 — 8' aperture Suzaku XIS source (NXB-subtracted) 
spectra and best simultaneous fit models to data from all three de- 
tectors (black data points and histogram in top panel), with resid- 
uals (bottom panel). Contributions from the ensemble of low mass 
X-ray binaries (if significant; LMXB, blue), from hot ISM thermal 
plasma components (red), and from the sum of GXB, CXB, and 
ICM background components (orange) are separately plotted. The 
NXB (purple), that is subtracted prior to spectral fitting, is also 
displayed. The top panel (a) shows the sum of the XIS0 and XIS3 
spectra, folded through a weighted response, the bottom panel (b) 
shows the XIS1 spectrum. 

cess absorption is strongly required (only) in the inner 4/ 

Ay 2 ~ 260 (26) in the 0 — 2' (2 — 4') region. The radial 
variation argues against the cause being incorrect model- 
ing of the contamination of the time-dependent Optical 
Blocking Filter, which we further confirmed by fitting 
the 0 — 2' spectra extracted from data reprocessed using 
a recently updated (2011 June) contaminant model. We 
have also applied the mixing model in Xspec (suzpsf) 
that aims to account for the cross-contamination of spec- 
tra due to the broad Suzaku point spread function (psf). 
These solutions tend to be unstable in how the ISM 
abundances are distributed; however, we note that the 
emission-averaged abundances over the entire (0 — 8') re- 
gion are robust and consistent with those derived from 
the above fits. 

We plot the Suzaku - derived abundance profiles of O, 
Ne, Mg, Si, S, and Fe in Figures 5abc, and abundance 
ratio profiles (with respect to Fe) in Figures 6ab. Since 
these profiles are not deprojected, nor corrected for the 
Suzaku point spread function, the true abundance gradi- 
ents are steeper. There is clearly a negative abundance 
gradient for Ne, Mg, Si, S, and Fe, and a hint of positive 
gradients in S/Fe and Si/Fe (Figure 6b). 

We plot the Suzaku - derived abundance ratios with re- 
spect to Fe for the entire galaxy (out to 8') in Figure 
7, estimated as an emission-measure- weighted average of 
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Table 4 

Suzaku Best-Fit ISM Parameters 


region 

xA l v 

N h 

kT 

0 

Ne 

Mg 

Si 

s 

Fe 

0 - 

2' 

1994/1624 

8-4to.6 

0-867/°°O3 

0 91+ 0 09 
-0.12 

2 09 + °- 19 
z - uy -o.n 

1 14 +o.05 
1 -- L - 0.04 

1 oi +0 ' 05 

1,U1 — 0.034 

1 20 + ° 15 
± - zu -o.io 

q 99+0.02 
u,yy -0.01 

2 - 

4' 

1578/1463 

io.otl;i 

0.890±°;°?° 

0 70+°- 17 
u - 0.17 

1 24+ 0 - 24 
0.18 

0 49+ 0 06 
u.^-0.07 

0 57+ 0 - 06 
u ‘°‘ -0.06 

0 75+ 016 
u - /o -0.16 

q 47+0.O8 

-0.04 

4- 

6' 

1529/1449 

2.2* 

0 931 +0.037 

u.»o j-_ 0 04 y 

0 9i+°- 42 
-0.31 

0 72 - *" 0,37 
u - ‘ -0.50 

0 35+ 0 - 20 
u ‘ t3D -0.15 

0 62 +0 ' 18 
u - oz — 0 .15 

0.56+« 

q 34+0.06 

U-o^—0 10 

6 - 

8' 

1054/1107 

2.2* 

0-87lt°;° 8 3 ? 

0-62±°-j 8 

0 046”*~°' 33 
U.U'TO — 0.046 

0 21”*” 015 
U,Z1 -0.16 

q 04+0.16 

u.o^_0 17 

0-38 ±S;|I 

n 1 c+0.09 
U-lo -0.02 

2 - 

6' 

2040/1946 

9 7+ 2 - b 
y./_ 2 7 

0-897i°;°i® 

0.84t°;« 

1 oi+°- 23 
1 - U1 — 0.16 

0 44+ 0 05 
u *^-0.07 

0 54+ 0 - 05 
u -°^-0 .06 

0-68±°;l® 

0 4i+ 0 06 

4- 

8' 

2019/2018 

2.2* 

0 923+ 0,025 
u.t/^o_ 0 062 

0.68t°;g 

0 31”*” 0 ' 17 
U '°‘ L — 0.22 

0 oc+0-14 

U-^o_Q 09 

n 44+0.17 
0.09 

o.43±g;5| 

n oo+0.02 

u,zz — 0.023 


Note. — Displayed for each spectral extraction region are minimized x 2 /dof for simultaneous fits to Suzaku XISO, XIS3, and XIS1 spectra, 
best-fit columns densities (in units of 10 20 cm -2 ), best-fit ISM temperatures (in keV), and best-fit O, Ne, Mg, Si, S, and Fe abundances referred 
to the solar standard of Asplund et al. (2009). Asterisks denote fixed parameters. 
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Figure 2. Same as Figure 1 for the 4 — 6 X aperture. 

the separate 0 — 2', 2 — 4', 4 — 6', and 4 — 8' spectra (con- 
sistent results are derived using the 0 — 2', 2 — 4', and 
4 — 8 ' regions, or the 0 — 2', 2 — 6', and 6 — 8' regions). 
In subsequent sections we use this as a proxy for the gas- 
density-averaged abundance. Analytic and semi-analytic 
estimates utilizing the /9-model fit to the surface density 
(see above), and a range of abundance gradients consis- 
tent with the data for various elements, show that the 
emission-measure approximation leads to overestimates 
of less than 25%. For comparison, the NGC 4472 pattern 
(renormalized to the same solar standard, but averaged 
over ~ 4.6i? e ) from Paper I is also shown. Although the 
average Fe abundance is 0.51 for NGC 4649 and 1.65 for 
NGC 4472, the Mg/Fe, Si/Fe, and S/Fe abundance ra- 
tios are consistent between the two galaxies. However, 
the agreement is poorer for the more problematic O/Fe 
and Ne/Fe ratios (O/Fe because of the importance of 
Galactic halo emission in the region of the O lines, Ne/Fe 
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Figure 3. Same as Figure 1 for the 2 — 4' aperture, 
because of blending of Ne with Fe L lines). 


3.2. XMM-Newton Spectral Analysis and Comparison 
with Suzaku 

Compared to Suzaku , we found the XMM-Newton 
spectral fits considerably more sensitive to the treat- 
ment of the (substantially more prominent) background 

- likely due to the presence of North Polar Spur emission 
(Randall, Sarazin, & Irwin 2004). In the end, we adopt 
what may be considered the most conservative approach 

- confining the fit to the 0.3-5. 5 keV band (MOS) or 
0.4-5. 5 keV band (PN) and, as detailed in the previous 
section, following the XMM-ESAS prescription that ex- 
clusively subtracts an estimate of the quiescent particle 
background. The XMM-ESAS generation of instrumen- 
tal backgrounds does not remove the fluorescent spec- 
tral line features; and, some residual soft proton con- 
tamination may also be present in the extracted spectra. 
Thus we include the following additional NXB compo- 
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Figure 4. Same as Figure 1 for the 0 — 2' aperture. 

nents in our EPIC spectral models: narrow, fixed-energy 
Gaussian emission lines corresponding to A1 Ka at 1.49 
keV (MOS and PN) and Si Ka at 1.75 keV (MOS only) 
to account for the former; and, a powerlaw with index 
fixed at 0.3 for the latter. As with the XIS spectral fits 
we include thermal components, modeled with the apec 
plasma code, to represent LHB, MWH, and ICM emis- 
sion, and a power law to account for the unresolved CXB. 
We again find that the addition of a narrow 1.24 keV line 
improved fits in the inner region, and this is included in 
0 — 2' and 2^4' EPIC, as well as RGS, fits. 

We display the best-fit model parameters for simulta- 
neous EPIC MOS1, MOS2, and PN 0 — 2' region spectral 
fits, and also for fits to the RGS spectra, for the com- 
bined observations (2001 and 2007) in Table 5. Spectra, 
best-fit models, and residuals are shown in Figure 8 for 
the former, and in Figure 9 for the latter. We compare 
0 — 2' Suzaku , 0 — 2' EPIC, and RGS abundance determi- 
nations in Figure 10, finding good agreement in general 
and excellent agreement with respect to the temperatures 
and Fe abundances. A discrepancy in the intrinsic ab- 
sorption persists, with a column density consistent with 
the Galactic value of 2.2 x 10 20 cm~ 2 (Dickey & Lockman 
1990) for the EPIC fits. We note that, other than the 
reduction in goodness-of-fit, fixing the column density in 
fits to Suzaku spectra at the lower Galactic value has 
little effect on the best-fit parameters - the largest ef- 
fect is a ~ 25% reduction in the 0 — 2' O abundance. 
Fixing the column density in the 0 — 2' EPIC spectra 
at the best-fit Suzaku value of 8.4 x 10 20 cm -2 does not 
impact the goodness-of-fit and, likewise, shifts best-fit 
abundance by ~ 20% or less. An even smaller effect is 
found if we restrict the XIS fits to energies >0.7 keV 
and fix the column density at the Galactic value. 





Figure 5. Fe and O ( top panel, (a)), Si and S ( middle panel, (b)), 
Mg and Ne ( bottom panel, (c)) abundance profiles. The points 
represented by circles are shifted to the right by 0.5 / . 

The results of fitting the EPIC spectra extracted from 
the 0—2', 2—4', and 4—6' annuli are displayed in Table 6. 
The higher background in the EPIC cameras did not al- 
low us to successfully fit the 6-8' annulus. Previous anal- 
ysis of XMM-Newton data (Randall, Sarazin, & Irwin 
2006; Nagino & Matsushita 2009; Ji et al. 2009) are gen- 
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Table 5 

XMM-Newton RGS and EPIC 0 — 2' Best-Fit ISM Thermal Parameters 


Spectra 

EPIC (2001+2007) 
RGS (2001+2007) 


X 2 + 

1756/1202 

822/753 


kT 

0-8621°;°“ 

0-8661°;“? 


O 

0 f>4+ 0 - 05 
U.O4_ 0 Q6 

n+0.15 

-0.15 


Ne 


l.CP 


1 . 66 ^ 


Mg 


I cn + 0.18 

±.oz_o i2 
o+0.34 
—0.33 


1 Q4+ 004 
1-04-0 05 

4+0.30 
*—0.27 


0.941 


Si 


1.25 

1.51 


+0.04 

0.05 

+0.60 

0.57 


1.41 


+ 0.11 

0.12 


Fe 


0.97l°;°i 

0.97+°“ 


Ni 


3.03 

2.93 


+0.14 

0.13 

+0.62 

-0.56 



u i i i i r 


T 


1 1 i 1 1 r 


«■> 2.5 

p-i 
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cd 
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CD 
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cd 
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Figure 6. Profiles of abundance ratios with respect to Fe for O 
and Mg ( top panel, (a)); and, Si and S ( bottom panel, (b). 

erally in agreement with our results although there are 
some significant discrepancies at some radii. As noted 
above (Section 3; see, also, Appendix A), comparison 
is complicated by the recent AtomDB update and, in 
fact, we find fair agreement in the inner regions if we 
shift to the older AtomDB 1.3.2 version and include a 
second, hotter ISM component (slightly favored in this 
case). Moreover, these studies employ simpler and less 
conservative background treatments that utilize subtrac- 


Figure 7. Abundance pattern comparison between NGC 4649 
(solid squares) and NGC 4472 (solid triangles). The NGC 4472 
values from Paper I have been rescaled to the Asplund et al. (2009) 
solar standard. 



CK5 ’ 1 2 ' ' 5 

Energy (keV) 


Figure 8. 0 — 2' 2001 XMM-Newton EPIC spectra (MOS1: black, 
MOS2: red, PN: green) and best-fit model with residuals (Table 
5). The spectra are fit over the 0.3-5. 5 keV band for the MOS and 
the 0.4-5. 5 keV band for the PN. 

tion of spectra extracted from “off-source” regions, while 
we explicitly include GXB and ICM components, allow- 
ing for their local variation and that of the instrumental 
background. For example, our central Fe and O abun- 
dances in the inner 2' agree well with previous work while 
we find lower Si abundance than (Nagino & Matsushita 
2009; Ji et al. 2009) in the same region. In the 2 — 4' and 
4 — 6' rings, where the background is more significant, 
we find higher abundances than these previous works. 

We plot the ratio of the EPIC-derived abundances to 
those of Suzaku in Figures llabc, and the EPIC abun- 
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Table 6 

XMM-Newton Annular Best-Fit ISM Parameters 


region 


kT 

0 

Ne 

Mg 

Si 

s 

Fe 

o 

1756/1202 

0-862/°;°“ 

0 64”*” 0,05 

1 r oTO. 18 
i.OZ—o 12 

i 04+0.04 
0.05 

1 or + 0.04 

I.zo_ 0 05 

1 41+o.u 
0.12 

0 97+ 0,01 
u - y ' -0.01 

2-4' 

1114/910 

0.966/°°?° 

l-14/an 

1 OO+0.71 

l.zO—o 24 

1 29” 1 ” 0,08 

1 46" 1 ” 008 

1 97+ 0.22 

-0.27 

1 09" 1 ” 0,01 
±.u»_o 03 

4-6' 

809/621 

0-960/°;°« 

0-50/°/° 

o.86/g;« 

0 3i+ 0 - 22 
u -° - 0.20 

0 70 +0 - 21 
u - ' u _0.12 

0 99+ 0,41 
u.yy_ 0 36 

0 45+ 0,12 


Note. — Displayed for each spectral extraction region are minimized x 2 /dof for simultaneous fits to XMM-Newton EPIC spectra. Shown are 
best-fit ISM temperatures (in keV), and best-fit O, Ne, Mg, Si, S, and Fe abundances referred to the solar standard of Asplund et al. (2009). The 
column density is fixed at the Galactic value, 2.2xl0 20 cm -2 . 



1 


1.5 


Energy (keV) 


Figure 9. RGS spectra (2007 Observation) and best-fit vapec 
thermal plasma model (Table 5), and residuals. RGS-1 ( top panel, 
(a)) and RGS-2 ( bottom panel, (b)) are broken out for clarity. First 
and second order spectra, included in the fitting, are shown. 


dance ratio profiles (with respect to Fe) in Figures 12ab. 
The latter show a dropoff in the 4 — 6' annulus with re- 
spect to the inner 4'. The EPIC S/Fe and Si/Fe (Figure 
lib) profiles echo those seen in the XIS data (Figure 6b) 
in indicating positive gradients. However, while XMM- 
Newton and Suzaku abundances show good agreement in 
the 0 — 2' and 4 — 6' regions, the EPIC abundances are 
more discrepant in the 2 — 4' annulus. We have inves- 
tigated this annulus in more detail by conducting joint 
EPIC /XIS fits, considering different bandpasses, and ex- 
amining fits with the absorption column density both free 
and fixed. In doing so we found that the EPIC spectra 
admit of two fits of comparable quality - one with Galac- 
tic absorption and relatively high abundances (Table 6) , 
one with elevated absorption and lower abundances that 
is in accord with the Suzkau fits (Table 4) . However the 
converse is not true: the XIS spectra are not consistent 
with the higher abundance model. 

Based on the above considerations, we adopt the abun- 



Figure 10. Abundance pattern comparison for separate fits to 
XMM-Newton EPIC 0—2' (solid squares) and RGS (solid triangles) 
and Suzaku XIS 0 — 2' (solid circles) spectra. 

dances that we infer from fitting the Suzaku spectra 
for comparison with simple chemical evolution models. 
These are evidently more robust, for the crude spatial 
binning that we employ, because of the lower background 
and sharper energy resolution of Suzaku (see Paper I). As 
a result we extend the measured abundance pattern to 
greater radii, encompassing a larger fraction of the opti- 
cal galaxy than was previously possible. A parallel anal- 
ysis that utilizes the XMM-Newton results in the same 
quantitative results and conclusions. 

4. MODELING THE NGC 4649 ABUNDANCE PATTERN 

Although the host of a weak AGN (Shurkin et al. 
2008; Ho et al. 2009; Dunn et al. 2010), NGC 4649 is not 
strongly disturbed optically (Ferrarese 2006) or in X-rays 
(Humphrey et al. 2008) and has an old, passively evolv- 
ing stellar population (Bregman et al. 2006). However, 
there are strong indications for complex time-dependent 
circulation flows in the NGC 4649 ISM. On one hand, 
the thermodynamic characteristics of the gas within 150 
pc and the flattened isophotes further out inferred from 
Chandra provide evidence for inflow from kpc scales all 
the way into the 3 x 10 9 M 0 (Humphrey et al. 2008) nu- 
clear black hole (Diehl & Statler 2007; Brighenti et al. 
2009). On the other, the fact that this material is rapidly 
cooling, yet has not appreciably accumulated in the core 
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arcminutes 

Figure 11. Ratio of the abundances derived from the XM Ad- 
Newton EPIC data to those of the Suzaku XIS detectors. We show 
Fe and O ( top panel, (a)), Si and S ( middle panel, (b)), Mg and Ne 
(bottom panel, (c)) ratio profiles. The points represented by circles 
are shifted to the right by 0.5 / . 

(in gaseous or stellar form), implies the existence of an 
outward transport mechanism that is probably intermit- 
tent (Brighenti et al. 2009). The impact of these flows on 



arcminutes 



arcminutes 

Figure 12. Profiles of abundance ratios with respect to Fe for O 
and Mg ( top panel, (a)); and, Si and S ( bottom panel, (b)) derived 
from the EPIC data. 


the integrated ISM abundances depends on the metallic- 
ity in the flow and its ultimate disposition - in particular 
whether it escapes the galaxy entirely and permanently. 

We utilize the abundance pattern derived from fits to 
the Suzaku spectra in constraining the enrichment of the 
hot ISM in NGC 4649. We adopt a global perspective 
in our interpretation of the NGC 4649 abundance pat- 
tern by seeking explanation of the hot ISM abundances 
measured over the entire optical galaxy via X-ray spec- 
troscopy as derived in previous sections, and defer de- 
tailed consideration of abundance, and possible abun- 
dance ratio, gradients. This is partially driven by the 
limitations of the Suzaku angular resolution, and the at- 
tendant degeneracies and uncertainties in spectral mod- 
eling, as well as the relative robustness of the galaxy- 
wide emission averaged abundances (and ratios). Our 
goal is to examine the overall sources and sinks of the 
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metals measured in this way, and attempt to quantify 
their relative contributions. In this way, we may identify 
signatures of the metal transport and enrichment mech- 
anisms discussed above, and better understand the “ins 
and outs” (Mathews 1990) of the gas flow in this elliptical 
galaxy. 

4.1. Conservation Equations and Steady-State Solution 

We proceed with a chemical evolution formalism along 
the lines of that in Paper I, with the notable addi- 
tions of explicitly allowing for inflow, and the inclusion 
of recently published SNIa yields (Maeda et al. 2010). 
The equations for the conservation of overall ISM mass, 
AfisM, and mass of the ith element, /ismMism, in a fully 
formed elliptical galaxy are as follows: 

— = M M r - A/out + M in , (1) 

dt 

and 

“ = ^SNIa+^MR - (l + ^out)A^out/lSM+A^in/h 1 ! 

(2) 

where /j SM is the mass fraction, Mmr and G^ R are 
the mass injection rates for the total gas mass and 
for the ith element from evolved stars, and G§ NIa the 
metal mass injection rate from SNIa for the ith ele- 
ment. G^ r = M M R/stars> where / s \ ars is the mass frac- 
tion of the ith element in these mass-losing stars, and 
G SNia = ^SNial/sNia’ where lV SN ia is the SNIa rate and 
j/gNia is the ith element SNIa mass yield. M 0 ut (Min) is 
the rate of flow out of (into) the system, b l out a bias factor 
allowing for an abundance offset for the outflowing gas 
(Paper I), and /* n the mass fraction of the ith element in 
the inflowing gas. Af ou t refers to the flow of gas beyond 
~ 6.5 R e: Mi n to the introduction of gas into this same 
volume by any means other than that associated with 
stellar mass loss. For a single-phase outflow one expects 
b l out < 0 as a result of negative abundance gradients in 
the flow, and 6* ut > 0 if gas of higher metallicity prefer- 
entially exits the galaxy due to buoyancy effects in hot 
SNIa ejecta (Tang & Wang 2010). For simplicity we as- 
sume that b l out = b out is identical for all elements (Paper 
I) , and that the abundances in the inflow are in the same 
proportions as in the stars, /* n = <&n/stars- 

In the present work we consider equations (1) and (2) 
in the steady state limit, in which case the source terms 
may be considered as averages over the density- weighted- 
mean residence time of the gas. Solving for the ISM 
abundances, 

fi _ /'stars ( ^ A Hn An ) A Q !J g \ [ a / 0 , 

/,su " (l + l>;„)(l+r,„) • (3) 

where r in = Mi n /M M R = M out /M M R ~ 1 and q = 
O.OO6676» S Nia/0MR Mq 1 . 6*mr is the specific mass loss 
rate in units of 2.4 x 10 -11 M Q yr _1 A R q (Paper I), 
and 0sNia is the specific SNIa rate in units of the es- 
timated present-day elliptical galaxy rate of 0.16 SNU 6 

6 SN per century per 10 10 Lbq 


(Cappellaro et al. 1997; Cappellaro, Evans, & Turatto 
1999; Dilday 2008). 

4.2. Source Terms 

ISM abundances depend on the ratio #SNia/0MR 
through the parameter q in equation (3). We generally 
assume $mr = 1 when quoting values of $sNia based on 
observationally inferred values of q. The corresponding 
stellar mass return rate (at the present time), derived in 
Paper I, is ~ 40 — 60% times higher than often assumed 
based on older estimates of turnoff and stellar rem- 
nant masses (Faber & Gallagher 1976; Mathews 1989; 
Ciotti, Pellegrini, & Renzini 1991). Mmr decreases with 
time approximately as f -1 - 3 (Paper I), so that the ef- 
fective value of 6*mr may be greater if some portion of 
the ISM was lost from stars at an earlier epoch - the 
average over the past 10 Gyr is 2.2 times the present 
value. At the same time, estimates of the number of 
SNIa produced per stellar mass formed (~ 0.0005 Mq ; 
e.g., Horiuchi and Beacom 2010) result in an aver- 
age SNIa rate three times our fiducial rate of 0.16 
SNU (though ~ half of these may be prompt explo- 
sions that are irrelevant for enriching the present-day 
ISM). Also, a recent estimate of the SNIa rate (Li et al. 
2011) is ~twice our fiducial value, and there is evi- 
dence of an enhanced rate in clusters (Mannucci et al. 
2008; Dilday 2010; Maoz, Sharon, & Gal-Yam 2010). 
Further complication arises due to the absence 
of a compelling theoretical explanation for rates 
as high as these (Ruiter, Belczynski, & Fryer 2009; 
Mennekens et al. 2010). In other words, dsNia = $mr = 
1 represents a reasonable benchmark but not a strict ex- 
pectation: the larger values estimated from some optical 
surveys lack a firm theoretical basis that would indicate 
the magnitude and evolution of the SNIa rate in any indi- 
vidual elliptical galaxy with its particular structure and 
stellar population. 

In Paper I, individual stellar abundances are de- 
termined with reference to a specific SNII yield set 
(Kobayashi et al. 2006) and the assumption that the 
SNIa enrichment of the stellar and ISM components 
are identical; the SNIa/SNII mix was then related to 
the stellar a/Fe ratio that was allowed to vary. Here 
we adopt the simpler assumption that stellar O, Ne, 
Mg, Al, Si, S, and Ar abundances are in solar propor- 
tions and overabundant with respect to Ca, Fe, and 
Ni by a fixed, constant factor given by the a/Fe ra- 
tio. We base this on the consistency we found in Pa- 
per I between the optically measured [a/Fe] s tars and 
those in our best-fit models to the hot ISM abundance 
pattern in NGC 4472. Moreover, there is no consen- 
sus as to which of the many available SNII yield sets (if 
any) is appropriate (Gibson, Loewenstein, & Mushotzky 
1997; Romano et al. 2010), especially given the uncer- 
tainties in, e.g., the initial mass function of the orig- 
inal stellar population and the possibility that SNII 
yields vary as the metallicity of the stellar popula- 
tion that provides the Type II progenitors accumu- 
lates. There is also the possibility that the primor- 
dial enrichment of the stars and the ongoing enrich- 
ment of the hot ISM are characterized by distinct 
yields from separate components of a bimodal popula- 
tions of SNIa (Mannucci, Della Valle, & Panagia 2006; 
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Maoz et al. 2011; Horiuchi & Beacom 2010) or, as may 
be the case with SNII, that SNIa yields evolve with the 
progenitor population - in their case due to effects of age 
and metallicity on the deflagration-to-detonation tran- 
sition density (Bravo et al. 2010; Krueger et al. 2010; 
Jackson et al. 2010). 

In this way, we focus on the period subsequent to the 
epoch of formation of the bulk of the stars, and only as- 
sume that an old stellar population with the optically 
estimated [a/Fe] stars has emerged during the unspeci- 
fied prior evolution. However, one must keep in mind 
that there is evidence that individual a-elements devi- 
ate from solar ratios in massive ellipticals, and do not 
vary in lockstep in general (Graves et al. 2007, Smith 
et al. 2009, and references therein) - in fact our group- 
ing of Ca with Fe above is based on just such evidence 
for that element (see, also, Pipino & Matteucci 2004; de 
Plaa et al. 2007). Following Humphrey & Buote (2006); 
Humphrey et al. (2008), we estimate the following global 
stellar values from measurements of optical line indices 
(Trager et al. 2000a, b; Thomas et al. 2005; Howell 2005; 
Sanchez-Blazquez et al. 2006): [Fe/H] s tar s = 0±0.2 and 
[a/Fe] s tars = 0.25 ± 0.05, and adopt [a/Fe] sta rs = 0.25. 
The only relevant supernova yields are those of SNIa that 
explode after long (>Gyr) evolutionary delay times. We 
consider a wide range of published SNIa yield sets, in- 
cluding those derived from the W7 (Nomoto et al. 1997) 
and more recent C-DEF (Maeda et al. 2010) pure defla- 
gration models, and those based on one-dimensional de- 
layed detonation (DD) models with a range of central and 
transition densities (Nomoto et al. 1997; Iwamoto et al. 
1999). For the first time, we consider the yields based on 
the two- dimensional DD models of Maeda et al. (2010), 
including their off-center DD model (O-DDT). 

Finally, in comparing our models with the NGC 4649 
ISM abundance pattern we consider only elements that 
have distinct and well-measured emission features in the 
X-ray spectra - O, Mg, Si, S, and Fe although we also 
plot Ne when displaying this comparison. 


4.3. Solutions Applied to NGC 4649 

We consider two classes of simplified solutions (with 
neither outflow nor inflow) to equation (3) for illustra- 
tive purposes. In the first the SNIa rate is set to 0 
( q = 0) so that ISM abundances are proportional to stel- 
lar abundances. Figure 13a includes a comparison of the 
hot ISM and stellar abundance patterns for solar and 
half-solar stellar Fe abundances ([Fe/H] sta rs = 0, -0.3). 
Also displayed is the best-fit model including SNIa with 
0SNia/0MR = 1: the [Fe/H] sta rs = -0.3 C-DEF model 
(where SNIa yields for relevant elements are lowest). No- 
tably, abundances are generally overpredicted - even for 
the [-Fe/H] stars = —0.3 model with no SNIa. 

We now turn to more general steady state solutions. 
Once the stellar abundance ratios and SNIa yields are 
specified, the ISM abundance ratios depend on a single 
parameter. Adopting Fe as a reference element, the ra- 
tios are determined by the following: 


/ism 

fFe 
J ISM 


/starsZ/stars ? PVsNla 

1 + PV SNIa 


(4) 


for each element i, where / s \ ars // s ^® rs is the correspond- 
ing stellar ratio (the solar ratio for Ca and Ni, and 


1.8 times the solar ratio otherwise; see above), and 
p = q{ 1 + Dn<Zin) _1 // S te r s- Essentially, this expresses 
the simple dependence of the ISM ratios on the ratio of 
SNIa-to-stellar direct injection, modified by inflow (for 
ri n > 0); a uniformly biased outflow does not affect the 
ratios. 

We determine the best-fit (minimum y 2 ) value of p 
to the observed NGC 4649 abundance ratios for each 
yield set. Figure 13b shows the corresponding abundance 
patterns for the more recently published yields compared 
to the observations. Many of the yield sets provide good 
fits to the abundance ratios (as does the W7 model); the 
exception being the Maeda et al. (2010) C-DDT model 
with its exceptionally high Si-to-Fe ratio. We do not 
find evidence for the oft-reported underabundance of O 
(see Paper I, and references therein) - in part because of 
the lower solar standard O abundances in Asplund et al. 
(2009) that we adopt. 

4.4. Discussion 

4.4.1. Astrophysical Interpretation of Successful Models 

How do the dimensionless model solutions that match 
the observed hot ISM abundances translate into astro- 
physical quantities of interest, and are these reasonable? 
We address that in this subsection by unfolding the sin- 
gle model parameter p into the implied range in the 
combination of inflow, outflow, and SNIa properties in 
NGC 4649. The value of p, determined by the observed 
abundance ratios, determines the relation between the 
inflow term and the SNIa-to-stellar-mass-return 

ratio term q (see above), and sets a minimum value of 
q {qi nin ) such that r in (j) in > 0: q m ; n = P f stars ' If 0ne 
then specifies one may separately derive rj n and 
< pi n for a given outflow bias parameter, & out . The min- 
imum value of q required for n n > 0 and </>i n > 0 is 
qmin = max[pf* t l rs , {l+b out )pf[£ M / {l+py sma )}. In what 
follows we set $mr = 1 when quoting values of $sNia (he., 
$SNia = 150g), and assume [Fe/H] sta rs = 0.0, unless oth- 
erwise specified. For illustrative purposes we adopt a 
composite model with equal contributions from C-DEF, 
O-DDT, and WDD3 SNIa (the three distinct well-fitting 
models in Figure 13b); the following discussion is qual- 
itatively unaltered, and the goodness-of-fit similar, for 
any of the individual well-fitting models. The abundance 
pattern inferred from optimizing the match to the data 
over the parameter p for this composite model is shown in 
Figure 14, where it is compared with the pure stellar and 
SNIa ratios, as well as the predicted ratio for $sNia = 1 
and no inflow {r\ n (fi n = 0). The theoretical pattern is 
extended to the additional elements C, N, Al, Cr, Mn, 
and Ni (stellar abundances of the former three are as- 
sumed equal to those of O-through-Ar, the latter three 
those of Ca and Fe; see Section 4.1.2). Figure 15a plots 
nn</>in as a function of 0sNia for this optimum composite 
model (solid line in Figure 14). The minimum SNIa rate 
(nn</>in — > 0), set by the lower value of [a/Fe] in the ISM 
compared to the stars, corresponds to 0SNia,min = 0.18 
(0.091 for [Fe/H] stars = -0.3); more generic.ally 0 SN ia > 
0.12 for all of the well-fitting individual SNIa yield sets 
(for [Fe/H] s tars = 0.0). The curves in Figure 15a may 
be expressed as (9 S Nia = 6 , SNia,min(nn<Zin + 1)- That is, 
very large inflow rates are required to reconcile the NGC 
4649 ISM abundances with 0sNia = 1- 


NGC 4649 Abundance Pattern 


13 



atomic number 



Figure 13. Top panel (a): Comparison of NGC 4649 abundance 
pattern (filled squares) and three simple (no outflow or inflow) 
model predictions: stars-only enrichment with [Fe/if] stars = 0 
and —0.3 (upper and lower solid lines), and the best- fit model with 
[Fe/H] stars = —0.3 and SNIa exploding at the standard rate (C- 
DEF yields, broken line). Bottom panel (b): Best- fit solutions to 
equation (4) to the NGC 4649 abundance ratio pattern (excluding 
Ne) for the following SNIa abundance sets: Maeda et al. (2010) C- 
DEF ( p = 1.25), C-DDT ( p = 1.72), and O-DDT ( p = 0.94) mod- 
els (solid, dotted, and short-dashed lines), Iwamoto et al. (1999) 
WDD 1 ( p = 0.94) and WDD3 (p = 0.63) models (long-dashed, 
dot-dashed lines; the WDD2 solution closely resembles that for 
the O-DDT yields). The solutions, which predict only relative 
abundances, are scaled to the observed ISM Fe abundance. 

In Figure 15b, rj n and (f>i n are separately broken out 
based on the measured ^half-solar global value of /jg e M 
derived from the Suzaku spectral fits. The inferred abun- 
dances in the accreted gas are < 0.28 times the (nearly 
solar) stellar value for $sNia < 1, reasonable for an in- 
tracluster or intragroup medium - and consistent with 
the outermost Suzaku annulus (6 — 8') abundance (Table 
3), while r; n > 2.1. The introduction of a biased out- 


flow that preferentially carries metals out of the galaxy 
naturally reduces the require inflow rates. For 6 out = 1 
(abundances in the outflow twice that in the stars), the 
inflow metallicity simply shifts to lower inflow rates (Fig- 
ure 15b) for fixed $sNia < 1 ri n — > (r; n — l)/2; more 
generally, r in -»■ (r in - & ou t)/(l + bout) (with a compen- 
sating increase in (f > since r- ln (j)i n is independent of 6 0 ut)- 
For [Fe/.ff] s tars = —0.3 the required level of inflow en- 
richment is higher for a given SNIa rate (Figure 15a) 
- r in is independent of [Fe/H] s t ars ; for [Fe/H] s tar s -t 
s f ars (f) > [(rj n (/>i n + l)([F’e/iF] stars /[F 1 e/lF] s tars) 

iKn 1 - 


The relative importance of dilution is expected to in- 
crease with radius as the timescale to replenish the ISM 
via stellar mass increases, with part of this dilution per- 
haps arising from an ab initio group-like orgin of the gas 
there. We do not account for this, but will address this 
with future multizone models. 



Figure 14. Best-fit model to the observed abundance pattern de- 
rived from Suzaku data (filled squares) assuming equal numbers of 
C-DEF, O-DDT, and WDD3 SNIa (solid line), compared with the 
ratios for pure stellar (dotted) and SNIa (long-dashed) enrichment 
and that for 0gNla = 1, rj n 0 i n = 0 (standard SNIa rate, no inflow; 
short-dashed line). 


4.4.2. NGC 447% an d Other Ellipticals (Briefly) Revisited 

As noted above (Figure 10), the abundance ratio pat- 
terns are generally consistent between NGC 4472 and 
NGC 4649. This suggests that the higher abundance 
of Fe in the former is not simply due to a higher SNIa 
rate - which would skew the abundance pattern - but, 
in the context of the models described above, may be 
explained by either more efficient removal of enriched 
gas or a higher average inflow rate of relatively unen- 
riched gas in the latter. In general, the level and pat- 
tern of abundances for other elliptical galaxies observed 
with Suzaku are near solar (Matsushita et al. 2007; 
Tawara et al. 2008; Komiyama et al. 2009; Sato et al. 
2009; Hayashi et al. 2009) - similar to what we report 
here for NGC 4649, indicating that our simple models 
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Figure 15. Top panel (a): n n </>i n vs. 0gNla for the best-fit model 
(solid line in Figure 14). Also shown, is the corresponding curve 
for [-Pe/l?] stars = —0.3 (dotted line). Bottom panel (b): </>; n vs. 
r; n implied by the observed /^ e M for the best-fit abundance ratio 
patterns assuming [Fe/H] stars = 0 and fe ou t = 0.0 (solid curve) or 
1.0 (broken curve). The corresponding 0gNla are shown by crosses; 
the first five values are the same for each curve. Abundance ratios 
are invariant along the lines in (a); abundances are invariant along 
the curves in (b) for the specified SNIa rates. 


tion would be helpful in this regard. A variation (or 
lack, thereof) in the magnitude of dilution from galaxy 
to galaxy constrains the nature of the inflow and its re- 
lation to elliptical galaxy evolution (see below). 

4.4.3. The Significance of Ni 

With the effects of SNIa diluted by inflow, it becomes 
more challenging to use ISM abundances in ellipticals to 
constrain SNIa yields. Even so, predictions for the pat- 
tern among the Fe group elements are sufficiently dis- 
tinctive that some signature remains, as we illustrate in 
Figure 16 for the models in Figure 13b. In particular, 
the Ni abundance (both its absolute value, and that rel- 
ative to Fe and other elements) is a diagnostic of both the 
level of SNIa enrichment and of the SNIa nucleosynthe- 
sis prescription (see also Paper I, PM11). Also plotted 
in Figure 16 are the Ni/Fe ratios for fits to RGS and XIS 
0 — 2' spectra where the Ni abundance is allowed to float. 
At the ~keV temperatures of elliptical galaxies, these are 
determined by fitting the shape of the spectrum in the 
region where various Fe L (and Ni L) features are gener- 
ally blended. As a result, there are concerns that these 
measurements may not be robust to assumptions and 
uncertainties in the plasma code atomic physics input 
implemented in spectral fitting packages such as Xspec. 
As such it will crucial to confirm the high Ni abundance 
estimated in this work, Paper I, Ji09, etc. with future 
high energy resolution X-ray imaging spectroscopic ob- 
servations, such as those to be conducted by Astro-H , and 
to improve the accuracy of thermal plasma modeling of 
Ni L. 
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may be widely applicable in accounting for the origins 
of ISM metals in elliptical galaxies (although with some- 
what lower abundances than typical of X-ray luminous 
ellipticals, NGC 4649 may be a somewhat extreme case). 
A more thorough analysis that compares and contrasts 
elliptical galaxy ISM abundances in the context of their 
environments and intrinsic properties is beyond the scope 
of this paper, but is currently being pursued. In partic- 
ular, a larger sample may help in evaluating which SNIa 
yield set is favored by the X-ray data (next section) . An 
abundance pattern that displayed a lower level of dilu- 


Figure 16. Same as Figure 13b, focusing on the Fe group ele- 
ments, with Ni/Fe abundance ratios from RGS (solid square) and 
XIS 0 — 2' (solid triangle) spectral fits. 


5. SUMMARY, IMPLICATIONS, AND FINAL REMARKS: 
RESULTS FROM X-RAY SPECTROSCOPY OF NGC 4649 
IN THE BROADER CONTEXT OF ELLIPTICAL 
GALAXY STUDIES 

The distinctive history of giant elliptical galaxies re- 
sults in the creation and maintenance of an extensive 
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ISM, dominated by hot gas, that contains a fossil record 
of that very history: evolution determines ecology, ecol- 
ogy enables archeology, archeology illuminates evolution. 
Based on estimates of supersolar stellar a-to-Fe ratios 
and SNIa rates in excess of > 0.1 SNU, one naively ex- 
pects an ISM abundance pattern in elliptical galaxies 
that is supersolar across the board, and increasingly so 
for elements more proficiently synthesized in SNIa. In 
actuality, measured ISM abundances are solar or sub- 
solar for all elements and do not deviate strongly from 
solar ratios, possibly excepting Ni. A simple reduction 
in the SNIa rate (or by assuming that SNIa ejecta fail to 
mix into the hot ISM; Brighenti & Mathews 2005) does 
not resolve this puzzle - the effective rate would need 
be reduced to very low values to explain the modest Fe 
abundance level, begging explanation for why the ratios 
do not then mirror the a-element enhanced pattern in 
the mass-losing stars. 

We bring these issues into sharper focus, and attempt 
to make progress on their resolution, through measure- 
ment of the abundance pattern in the elliptical galaxy 
NGC 4649 derived from deep Suzaku and XMM-Newton 
observations. The measurement of abundances in the 
hot ISM of elliptical galaxies remains problematic, de- 
spite the high quality of the spectra extracted from these 
data due to effects of angular resolution in the case of 
Suzaku, background systematics in the case of XMM- 
Newton, and lingering atomic physics uncertainties (see 
Appendix A). Nevertheless, by demanding a degree of 
cross-mission consistency, and focusing most on conclu- 
sions that rely on relatively model-independent abun- 
dance ratios integrated over the galaxy, we find that ro- 
bust constraints on elliptical galaxy evolution can be in- 
ferred. Towards this end, we compare the galaxy-wide 
average ISM abundance pattern with one-zone steady 
state solutions to the equations of chemical evolution 
that depend on a single parameter that characterizes the 
relative contributions of SNIa, stellar mass loss, and in- 
flow to the ISM metal inventory. We vary this parameter 
to fit the data, utilizing the most recent solar standard 
abundance scale (Asplund et al. 2009) and SNIa yields 
(Maeda et al. 2010), and reproduce the observed abun- 
dance pattern in NGC 4649, thus lending some insights 
into the contributions of various sources of ISM metal 
enrichment. 

Our approach has its drawbacks and limitations, and 
should be considered provisional as we expand the data 
analysis sample and refine our formalism. By adopting 
the steady-state limit of equations (1) and (2), we focus 
on the relative contributions of stellar mass loss, SNIa, 
and inflow/outflow (in the broad sense in which they are 
defined here), but gloss over the fact that these are not 
strictly co-eval and that the dilution factor refered to as 
“inflow” may in part correspond to a pre-exisiting group 
environment that continues to be reflected in gas phase 
abundances at large galactic radii. Encapsulating the 
relative impact of stellar mass return, dilution, and SNIa 
on the pattern of abundance ratios in a single parameter 
has the virtue of facilitating accurate model/data com- 
parison, and emphasizes their interplay and degeneracy 
- e.g., the challenge in distinguishing the combination of 
low SNIa rate and low dilution from that of high SNIa 
rate and high dilution. However, this complicates the 
interpretational framework somewhat while, at the same 


time, limiting it by leaving these unbroken degeneracies 
in place. We are confident that these degeneracies may 
be broken as the sample, and our confidence in absolute 
abundance measurements, grow. 

We find that a key ingredient in models that success- 
fully reproduce the observed pattern is dilution via in- 
flow of subsolar metallicity extragalactic gas with an (as- 
sumed) abundance pattern similar to that in the stars at 
an average rate comparable to, or greater than, that of 
stellar mass loss. As a result the overall ISM metallicity 
is reduced while, at the same time, the effectiveness of 
Fe group enhancement from direct SNIa injection is di- 
luted. Although SNIa rates of 0.1 SNU or more are not 
required, we now find that they may be accommodated 
(and, in fact, that specific rates < 0.02 SNU are excluded 
by limits on the ISM [a/Fe] ratio relative to that in 
the stars) provided that a significant fraction of the ISM 
is of external origin. Astrophysically reasonable magni- 
tudes of SNIa enrichment and external dilution may play 
off against one another, as shown in Figure 15, to pro- 
duce subsolar abundances in roughly solar proportions 
as observed. In a similar vein, Brighenti et al. (2009) 
modeled the two-dimensional gasdynamics of NGC 4649 
and found that inflow of circumgalactic gas with metal- 
licity comparable to the mean stellar value could di- 
lute the effects of SNIa enrichment (in their models, 
^SNia/^MR = 0.84 using our scale) on the Fe abundance 
by a factor of ~ 3. We cannot precisely identify the 
origin and means of delivery of the extragalactic mate- 
rial. Inflow may be quasi-steady, or intermittent follow- 
ing the episodes of outflow that must occur to prevent the 
gas from accumulating to levels beyond what is observed 
(Colavitti, Pipino, & Matteucci 2009). 

The reservoir from which gas accretes may be an intra- 
cluster, intragroup, or filamentary intergalactic medium; 
and may have been previously ejected from the same 
galaxy (or other galaxies). Alternatively, the extragalac- 
tic material may originate in discrete instances of galaxy 
merging with gas-rich systems. Some combination of 
these mechanisms would seem most likely, with their rel- 
ative importance depending on galaxy history and en- 
vironment (La Barbera et al. 2011). The level and pat- 
tern of ISM metal enrichment will reflect this and, with 
detailed modeling and additional observational analysis, 
may prove an important diagnostic of the nature of post- 
formation elliptical galaxy interactions. 

Thus, the X-ray abundance pattern confirms the com- 
plex flow and enrichment history of NGC 4649 pro- 
posed in Brighenti et al. (2009), and reflects the dy- 
namic, continually evolving nature of elliptical galax- 
ies implied by optical observations (Faber et al. 2007). 
The earliest epochs of elliptical galaxy formation are 
characterized by the interplay of major mergers, in- 
flow of primordial gas, vigorous star formation, and 
powerful outflows (Colavitti, Pipino, & Matteucci 2009; 
Mannucci et al. 2010; Arrigoni et al. 2010). Even af- 
ter the quenching of star formation has halted the 
main formation epoch, ellipticals do not evolve purely 
in a passive manner in isolation (Thomas et al. 2010). 
Evidence of mergers and other interactions, indicate 
that ellipticals continue to grow up to the present day 
(Trager et al. 2000b; van Dokkum et al. 2008; Tal et aL 
2009), as required to explain the observed increase 
in the size of ellipticals and the continual build 
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up of the red sequence (Naab, Johansson, & Ostriker 
2009; van der Wei et al. 2009; Shankar et al. 2010; 
Robaina et al. 2010; Trujillo, Ferreras, & de La Rosa 
2011; Damjanov et al. 2011; Cassata et al. 2011). In ad- 
dition, there are signs of the presence of modest amounts 
of more recently formed stars - often in ellipticals with 
morphological signs of interaction (Trager et al. 2000b; 
Kaviraj et al. 2008, 2009; Sanchez- Blazquez et al. 2009; 
Kaviraj 2010; Vega et al. 2010). However, evolution sub- 
sequent to the primary star forming epoch must proceed 
in such a way as to be consistent with the old and pas- 
sively evolving nature of elliptical galaxy stellar popula- 
tions inferred from their individual properties and scaling 
relations (Pipino & Matteucci 2006, 2008). Attempts to 
resolve this apparent paradox, alluded to in the intro- 
duction, often appeal to the process of “dry” merging 
(Tal et al. 2009; Kaviraj 2010; Cooper et al. 2012) that 
involves very little gas - and hence new star formation 
(as contrasted with the z > 2 “wet” mergers where stars 
are efficiently formed). The observed scaling relations 
and predominance of old stellar populations are thus 
preserved. However it is not clear that mergers of this 
type occur in sufficient numbers (Hopkins et al. 2010); 
and, most galaxies interacting with ellipticals are not, in 
fact, gas-free (De Propris et al. 2010; Serra & Oosterloo 
2010 ). 

Massive elliptical galaxies are themselves not gas-poor, 
and the observed properties of the hot ISM are often ne- 
glected in these considerations. We have shown that the 
chemical structure of the hot ISM in ellipticals not only 
confirms a continual environmental impact, but demon- 
strates that the accreted gas is hot or is heated as it is 
introduced and interacts with pre-existing hot gas. This 
implies a reduction of the star formation efficiency of any 
cold gas that may be newly introduced either in the form 
of mergers with small gas-rich galaxies or via a smoother 
accretion process (Keres et al. 2009). 

The entire one-dimensional “humidity-based” classifi- 
cation of mergers 7 implicitly assumes a one-to-one cor- 
respondence between gas content, dissipation, and star 
formation that breaks down in giant ellipticals known 
to be filled with hot gas. An encounter involving such 
a galaxy is never gas-poor, and may involve substan- 
tial dissipation if the companion galaxy is also gas-rich. 
However, these mergers may prove to be a hostile setting 
for new star formation, and are effectively “dry.” Alter- 
natively, the late-time accretion may be dominated by 
hot mode accretion (van de Voort, et al. 2011). Under- 
standing elliptical galaxy evolution requires tracking the 
three “phases” of baryonic matter - stars, (cold) star- 
forming gas, and (hot) inert gas - and how their mutual 
exchange of mass, metals, and energy as a function of 
age, environment, and type of merger proceeds as the 
galaxy evolves subject to internal and external influences 
(De Lucia et al. 2010; Lu et al. 2011). 

The rapid build-up of the bulk of the stars that 
compose present-day giant elliptical galaxies was 
quenched by some process that maintained a galaxy 


ecology free of significant amounts of cold gas. Suggesed 
mechanisms, that also may account for the observe 
dichotomy in the characteristics of galaxies, include 
quenching via gravitational heating (Dekel & Birnboim 
2008; Birnboim & Dekel 2011; Nipoti & Binney 2007; 
Johansson, Naab, & Ostriker 2009; Cattaneo et al. 
2009), and via AGN feedback (Bell et al. 2012 and 
references therein). We see that ISM abudances are a 
potential diagnostic of the continuing operation of star 
formation quenching. 

Our conclusions are driven by the robust analysis re- 
sults indicating that ISM abundances are solar or less 
with ratios that lie between the a-element enriched pat- 
tern expected from pure enrichment from stellar mass 
loss and the O/Ne/Mg-poor pattern expected from the 
combination of stellar mass loss and direct SNIa injec- 
tion at $sNia ~ 1. As NGC 4649 is a member of the 
Virgo Cluster and displays a compact X-ray morphol- 
ogy, the evolution of its hot ISM (e.g., stripping of the 
outer regions) may be affected by the interaction of the 
galaxy (or a subgroup to which it belongs) with the ICM. 
Nevertheless, given the universality of the departure of 
observed elliptical galaxy hot ISM abundances from the 
naive expectations described above, we suggest that the 
following main conclusions we draw for NGC 4649 may 
be generalized: 

• ISM abundance patterns indicate that metal en- 
richment from stellar mass loss and from direct 
injection of SNIa, and dilution from infall are all 
significant when averaged over the residence time 
of the ISM. SNIa injection prevents the ISM from 
reflecting the supersolar stellar [a/Fe] ratio, and 
infall dilutes the overall ISM metallicity. 

• An amount of low metallicity gas exceeding that 
originating in stellar mass loss rate must be intro- 
duced into the galaxy if one is to reconcile the ISM 
Fe abundance with the standard specific SNIa rate 
of 0.16 SNU. 

• We find further support to the picture of ellipticals 
as “open ecosystems” that continually grow, and 
exchange mass and metals with the intergalactic 
environment and/or other galaxies. 

• The gas that encroaches into or is introduced into 
the ISM must already be hot, or must be heated by 
mixing with the hot gas and/or some other internal 
mechanism in order to explain why these processes 
are accompanied by little or no star formation. 


We are grateful to Adam Foster, and Randall Smith 
for their input regarding atomic physics issues. This 
work could not have been completed without support 
from NASA ADAP , and Suzaku and XMM-Newton GO 
programs, and was improved by feedback from a con- 
structive referee’s report. 
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APPENDIX 

EVERYTHING YOU KNOW IS WRONG: EFFECTS OF UPDATING THE ATOMIC DATABASE 

We initiated analysis of NGC 4649 utilizing the AtomDB vl.3.2 data base as implemented in Xspec version 12.6. 
Subsequently, with the availability of the updated AtomDB v2.0.1 it became evident that the previous best fit models 
were no longer statistically acceptable due to severe differences in Fe L-shell strength predictions. This is primarily due 
to differences in the calculated ionization balance of Fe (Adam Foster, private communication). In this appendix we 
conduct simple illustrative numerical experiments to explore the effect of using the “wrong” database to fit spectra, at 
CCD energy resolution, of thermal plasmas with kT ~ 1 keV. We find significant discrepancies for kT < 1.3 keV that - 
to the extent that the update represents a better description of nature - casts doubt on the accuracy of previous X-ray 
results for elliptical galaxies and other ^keV astrophysical plasmas. The resulting artifacts include mis-estimates of 
the plasma temperature and metal abundances, and false statistical preference for two temperatures. 

First we construct simulated spectra, using the NGC 4649 0 — 2' background and response Hies. These spectra 
are normalized to yield 10000 source counts, and are modeled as a solar abundance (Asplund et al. 2009), single- 
temperature thermal plasma according to AtomDB v2.0.1. These are then fit with (one- and two-temperature) models 
that instead utilize AtomDB vl.3.2, and 90% confidence regions derived by cycling through 1000 such realizations. As 
a control, the identical experiment is conducted by repeating with AtomDB v2.0.1 fits. Temperatures of 0. 3-2.0 keV 
are considered. 

The ratios, with respect to the input temperature kTi Ut 2.0.I) of the best-fit AtomDB vl.3.2 temperature kT out ^\. 3.2, 
are shown in Figure 17a; kT outt 1.3.2 versus the best-fit control temperature fcT outi 2.0.1 in Figure 17b. Figure 18 shows 
the ratio of the best-fit Fe abundance in the AtomDB vl.3.2 fits to those in the control fits. Other elements show 
similar behavior, with larger errorbars. Clearly, temperatures and abundances are underestimated relative to their 
“true” values, with the effect most prominent between 0.5 and 1 keV. 

In addition, the goodness-of-fit statistic using the “incorrect” plasma model is significantly worse relative to the 
control fits for kT < 1.3 keV, as shown in Figure 19a. Given that the resulting reduced-y 2 (x 2 / v ) generally significantly 
exceeds one, standard procedure is to add an additional component. We conduct an experiment along these lines as 
well for simulation input temperatures of 0.8, 0.9, 1.0, 1.1 keV and find that, indeed, x 2 /^ generally declines (except 
for fcTjn^.o.i = 1-1) and attains acceptable magnitudes ~ 1.1 when an additional hotter (~ 1.2 — 1.6) temperature 
component is considered. This also converts what had been an abundance underestimate for one-temperature models 
using the “wrong” plasma model into an abundance overestimate , as demonstrated for Fe in Figure 19b. 



2.0.1 (keV) 
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Figure 17. Left panel (a): Ratio of the best-fit temperature using AtomDB vl.3.2, to the AtomDB v2.0.1 simulation input temperature, 
versus simulation input temperature. Right panel (b): Best-fit temperature using AtomDB vl.3.2. versus best-fit temperature using 
AtomDB v2.0.1. 

Based on these simple, illustrative numerical experiences we may conclude that deficiencies in the earlier versions of 
AtomDB that includes the apec thermal plasma model output spectra, likely resulted in temperature underestimates 
by as much as 15%. Metal abundance estimates were underestimated for single-temperature models, and overestimated 
for two-temperature models - where the latter would have been artificially preferred based on goodness-of-fit statistics. 
The effect, while not monotonic with temperature, is most significant for “true” temperatures ~0.5-l keV. 

Of course, while AtomDB v2.0.1 is surely more accurate than its predecessors various, instances of incompleteness 
and approximation necessarily remain, 8 begging the question of the accuracy of model parameters obtained in current 
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Figure 18. Ratio of the best-fit Fe abundance using AtomDB vl.3.2, to the best-fit that using AtomDB v2.0.1, versus simulation input 
temperature. 
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Figure 19. Left panel (a): Reduced-y 2 values for fits to simulated AtomDB v2.0.1 spectra with l:'71 rn ..2.a. 1 for fits using (the “correct”) 
AtomDB v2.0.1 single-temperature models (filled triangles), AtomDB vl.3.2 single-temperature models (crosses), and AtomDB vl.3.2 
two-temperature models (filled squares). Right panel (b): Same as Figure 19a for the Fe abundance, showing single- (crosses) and two- 
temperature (filled squares) models. 

spectral fits. 

When models of astrophysical plasma spectra are formally unacceptable, this may be attributed to complexities in 
the physical conditions or to inaccuracies in the plasma model. Distinguishing these explanations is facilitated by 
improvements in the quality of atomic databases and X-ray spectra. At the same time, confrontation between the 
models and data can point the way to improvements in the accuracy of both our understanding of the plasma physical 
conditions and the key remaining deficiencies in the atomic database and its application. 

7 in addition to “wet” and “dry” , mergers have also been referred to as “moist” (Sanchez-Blazquez et al. 2009) and “damp” (Forbes et al. 
2007) 

8 http://www.atomdb.org/ 
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